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A new form of nanoporous material, metal intercalated covalent organic framework (MCOF) is proposed 
and its energy storage property revealed. Employing density functional and thermodynamical analysis, we 
find that stable, chemically active, porous materials could form by stacking covalent organic framework 
(COF) layers with metals as a gluing agent. Metal acts as active sites, while its aggregation is suppressed by a 
binding energy significantly larger than the corresponding cohesive energy of bulk metals. Two important 
parameters, metal binding and metal-metal separation, are tuned by selecting suitable building blocks and 
linkers when constructing COF layers. Systematic searches among a variety of elements and organic 
molecules identify Ca-intercalated COF with diphenylethyne units as optimal material for H 2 storage, 
reaching a striking gravimetric density ~ 5 wt% at near-ambient conditions (300 K, 20 bar), in comparison 
to < 0.1 wt% for bare COF-1 under the same condition. 

Finding low cost, safe, and efficient energy storage materials has been a major challenge to develop renewable- 
energy based economy. Hydrogen is considered as an appealing energy carrier alternative to fossil fuels in 
many applications, because it is abundant, energy intensive, and pollution-free 1 " 3 . U.S. Department of 
Energy (DOE) set a target for practical H storage materials to have a gravimetric capacity of 5.5 wt% by 20 17 4 . 
The challenge is to dramatically increase storage capacity and reversibility at ambient conditions 5 . Carbon-based 
nanostructures and porous materials including fullerene, nanotube, and carbyne networks, have been proposed 
as potential H 2 storage media, thanks to their high surface-to-weight ratios, good reversibility and fast kinetics 6 " 10 . 
However, most of them failed at room temperature and ambient pressure due to improper binding to H 2 . 

To reach high and retrievable H storage at ambient conditions, a light three-dimensional (3D) material with H 2 
binding energies in the range of 0.2-0.6 eV/H 2 is required 11,12 . This implies porous materials must have dense 
active sites for H 2 binding, to which two approaches have been adopted. One is to substitute carbon by other 
elements such as boron/nitrogen 1113 ; however, it shows little effect and is difficult to apply. The other is to dope 
light metals including alkali, alkali earth or transition metals, into porous materials 1415 . Improving metal binding 
by using nanostructure curvature 1617 or incorporation of sp-carbons 1819 has been also explored. It was found 
alkali metals have too weak interaction with the substrate and H 2 20 ; whereas transition metals tend to cluster on 
material surface, reducing substantially H 2 storage capacity 21 . Upon compromise, calcium is recognized as a 
superior coating element thanks to its delicate balance of binding strengths to base materials and H 2 18,22 " 24 . Ca- 
coated fullerene, Ca 32 C 60 , yields strong local electric fields on C 60 and reaches a H 2 uptake up to ~ 8.4 wt% 22 . Ca- 
doped carbyne networks maintain a small distance ~ 10 A between metal adsorbates and a high H 2 density of ~ 8 
wt% 18 . However, both materials have their own flaws: their 3D assembly without affecting H 2 density is either not 
achievable or not experimentally demonstrated; and it is difficult to control the position and stability of dopant Ca 
atoms 18 ' 22 . Similar problems emerge in Ca-doped nanotube, graphite and graphene 14 ' 23 ' 24 . 

Covalent organic framework (COF) has been recently synthesized as an important extension of carbon porous 
materials 25 ' 26 . A variety of experimental studies 27 ' 28 show that small planar organic molecules form two- 
dimensional (2D) or 3D supramolecular architectures including COFs, through self-assembly of end linker 
groups, as schematically shown in Figure 1. Due to their adjustable porous structure and unparalleled strength, 
COFs are very promising for H 2 storage. At low temperature 77 K and high pressure 100 bar, 2D or 3D COFs can 
currently store 3-7 wt% H 2 27 . Stunning H 2 uptake of ~ 20 wt% at 77 K and 100 bar has been predicted in 
theoretical simulations for COF- 108 and COF- 102-3; but it deteriorates quickly to a value below 1-2 wt% at near 
ambient condition, for instance, at 300 K and 20 bar 29 " 31 . Similar trends were observed for COF-1, a simple layered 
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Figure 1 | Formation of COFs. (a) Schematic diagram of COF building block composed by two ending linkers and an organic unit, (b) An example of real 
building blocks: benzene diboronic acid (BDBA), (c) 2D COF formed by BDBA self-assembly, named COF-1. Crosses mark positions tested for metal 
binding. 



COF with the phenyl group as organic unit and boronic acid as the 
linker (see Fig. 1), shows a H 2 density of 1.28 wt% at 77 K but it 
decreases to ~ 0.08 wt% at 298 K and 20 bar in experiment 27,32 . 
Suitable metal doping also improves H 2 storage in COFs 33 " 38 . 
Grand canonical Monte Carlo simulations reveal that H storage 
capacity in Li-coated COF-202 can reach 4.39 wt% at 298 K and 
100 bar, which is three times higher than in undoped COF-202 34 . 
However, with the increase of Li concentration, Li atoms form clus- 
ters and H storage density reduces 35 . Using density functional theory 
(DFT) calculations, Wu et al. found that both Li and Ca forms clus- 
ters on COF- 10, which are 0.17 and 0.20 eV/atom more stable than 
scattered atoms, respectively. This reduces storage density from 3 to 
1.5 H 2 per metal atom. Ca dopants provide a stronger bonding to H 2 
(energy ~ 0.21 eV) compared with Li dopants (0.12 eV), due to more 
charge transfer from Ca to the substrate 35 . As a remedy, two-step 
doping, namely, metal decoration on B-substituted COF fragments, 
was also proposed 36 . On all these successful 3D reference materials, 
H 2 density is below 1-2 wt% at near ambient conditions 30 ' 31 . New 
methods and concepts for further improvement without metal 
aggregation are key to develop next generation H storage media. 

Here we present a radically different approach to improve H stor- 
age in COFs. We explore the feasibility of using a new class of 3D 
network — metal-intercalated COF (MCOF), for near ambient H 2 
storage. Metal binding in such new class of 3D materials (see the 
inset of Fig. 2a) has a binding energy significantly larger than that in 
bulk metals, thus intrinsically blocking metal segregation. Quantum 
mechanical calculations and thermodynamical analysis reveal that 
optimized MCOF material stores 4.87 wt% of H 2 at 300 K and 20 bar, 
greatly enhanced from the value of < 0.1 wt% for the bare COF-1 in 
the same condition 32 . 

The paper is organized as follows. After this introduction section, 
we first expound the idea and important considerations in construct- 
ing MCOF in the "Results" section. The feasibility of various light 
metal elements intercalating in model COF layers, COF-1, is tested. 
Ca is identified to be the most suitable gluing metal element. To 
further improve H storage performance, we start COF structural 
optimization through the selection of organic building blocks and 
linkers. In optimizing organic blocks, small organic molecules as well 
as polycyclic aromatic hydrocarbons (PAHs) are considered and 
compared to the reference molecule, benzene. We also investigate 
the effects of linker groups, where we find -B(OH) 2 linkers preserve 
H 2 binding properties of doped organic units. Ca intercalated COF 
layers with diphenylethyne organic units (CaCOF-d) are identified as 
an outstanding H storage material. After discussing H 2 storage prop- 
erties in 3D CaCOF-d layers at saturation and near ambient condi- 
tions, we give a brief discussion on how MCOF can be synthesized 
experimentally, and present a summary on our major findings. 



Results 

It seems natural to combine light metals and extended COF layers for 
near ambient H 2 storage. The two are both light and stable, and can 
potentially form new framework with dense active sites. In early 
experiments the newly synthesized COF materials are found to inter- 
act weakly with H 2 molecules, leading to low storage density at 
ambient conditions. In the meantime, various theoretical efforts 
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Figure 2 | Partial density of states (PDOS) projected onto the COF 
component (black lines) and the c/-orbitals of Ca (red lines) in (a) CaCOF- 
1 and (b) CaCOF-1 after H 2 adsorption. The PDOS projected onto H 2 
blue line, multiplied by 20 times) is also shown in (b). The inset in (a) 
shows the side view of CaCOF-1, and the inset in (b) shows the charge 
density associated with the states in the shaded region. 
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Table 1 | Metal binding energies (Em) and interlayer distances (L) 
for 3 D stacked MCOFs. The cohesive energies (E c ) for bulk metals 
are also listed for comparison 
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Bulk metal 
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MA) 


E M (eV) 


Ec (eV) 


Li 


3.7 


2.289 


3.9 


2.231 


1.605 


Na 


4.6 


1.408 


4.7 


1.490 


1.048 


Mg 


4.3 


0.057 


3.7 


0.202 


1.296 


K 


5.4 


1.596 


5.4 


1.708 


0.798 


Ca 


4.4 


2.830 


4.4 


2.786 


1.877 


Ti 


3.6 


5.058 


3.6 


4.993 


5.207 



demonstrated that metal doping cannot effectively improve H 2 stor- 
age due to metal clustering problem on single COF layers. We start a 
new approach to explore that whether metal atoms could be inter- 
calated between COF layers to form stable 3D stacks without metal 
segregation. The choice of workable metal elements and COF mate- 
rials should be tested to satisfy the requirements of being stable, light, 
and active for H 2 storage. 

Constructing metal intercalated COF for H 2 storage. We first take 
one of the simplest COF materials, COF-1, whose organic unit is a 
benzene ring and linker boronic acid, as our representative model 
substrate and investigate its interaction with various metal elements. 
A range of light metals including Li, Na, Mg, K, Ca, and Ti are 
considered, since they are most commonly adopted to decorate 
porous nanomaterials serving as active sites for H 2 binding. We 
first test the adsorption site of metal atoms on a single layer COF- 
1 . We find metal atoms prefer to adsorb onto the center of the organic 
unit, a phenyl ring, rather than on top of B, O, or B 3 0 3 rings. For 
instance, the binding energy of Ca (E Ca ) on all sites falls in the range 
of 0.001-0.008 eV, except on the center of phenyl ring, which is 0.77 
eV. This site preference is a result of Ca interacting with the n orbital 
of phenyl ring, which is absent from the electron rich boroxine ring. 
As shown later, the energy for Ca binding to the phenyl group in COF 
is an order of magnitude larger than that on isolated benzene (0.07 
eV), because the n orbitals are extended onto neighboring B atoms 
connected to phenyl rings in COF and Ca binding is thus enhanced 
by delocalized n electrons, as clearly displayed in Fig. 2b inset. 

Most importantly, the Ca binding energy increases to 2.83 eV 
when Ca and COF-1 layers assemble into multilayered metal inter- 
calation structure, the MCOF. The interlayer distance increases from 
3.35 A to 4.40 A upon metal intercalation and the stacking manner 
might also be changed. Ca binding on boroxine rings is ~ 1 eV less 
stable than on phenyl rings in intercalated 3D COF layers. We expect 
metal intercalation takes place simultaneously or with a small barrier 
through the pores of a diameter of ~ 10 A present in COF-1. The Ca 
binding energy is 3.37 eV if Van der Waals density functional 
(VDW-DF) 39 ' 40 is used to incorporate van der Waals (VDW) inter- 
actions; the additional 0.54 eV is attributed to the VDW interactions 
between adjacent COF layers. For comparison, the cohesive energy in 
bulk Ca is 1.88 eV/Ca in our calculation and the experimental value is 
1.86 eV/Ca 40 , much smaller than Ca binding energy in COFs. The 
energy for a second Ca binding around the Ca on the phenyl ring 
center is 1.81 eV with a Ca-Ca distance of 3.28 A; the energy is 
smaller than the binding energy of first Ca and the cohesive energy 
of Ca bulk, showing straightforwardly that Ca clustering is not 
favored. Vibrational analysis also confirms that no major negative 
frequencies occur in the vibrational spectrum of CaCOF-1. 
Therefore, Ca and COF-1 forms stable 3D materials (referred to as 
CaCOF-1) via metal intercalation. 

Some other metals also work in the same way as Ca to form 
MCOFs. In Table 1, we list the calculated parameters including the 
interlayer distance (L) and metal binding energy (£ M ) for metal 
intercalation COF-1. Metal insertion into the COF layers signifi- 
cantly expands the equilibrium interlayer distances, from L = 3.4 A 
for bare COFs 25 to L = 5.4 A for KCOF-1. The binding energies of 
metals to COFs vary from metal to metal, with E M for Mg being 
smallest (< 0.2 eV) and E M for Ti being largest (> 5.0 eV) in our 
calculations. This is in accord with the rather inert (active) nature of 
Mg (Ti) metals. However, both cases have the metal binding energy, 
£ M , smaller than the cohesive energy in bulk metals, £ c , as is the case 
for most metal doped nanostructures previously studied 22 . Sur- 
prisingly, on the other hand, Li, Na, K and Ca metals have E M 
significantly larger than the corresponding E c (see Table 1). This 
behavior renders metal segregation in such MCOFs thermodynami- 
cally unstable. Consequently, 3D network formed in this way is a new 
class of stable materials without any metal clustering problems. 



The new MCOF material can be approximately viewed as reduced 
intercalation compounds with metal dopants in cation states. Our 
Bader charge analysis shows that there are about 1.42 electrons trans- 
ferred from Ca to COF-1 during CaCOF-1 formation, and 0.84 elec- 
trons transferred from Li to COF-1. This new class of materials 
contains porous structures with stable, active metal sites, which pro- 
vide extraordinary advantages for molecular adsorption and storage, 
in particular, for H 2 storage. In this aspect, MCOF as the new class of 
3D material is significantly different from earlier proposals including 
those employing Ca 32 C 60 and other metal-doped materials, in that: 

i) It is an extended three-dimensional framework that can be 
readily exploited in practical applications; while Ca 32 C 60 is finite 
whose 3D assembly is not feasible without significantly reducing 
H 2 density. 

ii) MCOF is thermodynamically stable with metal binding ener- 
gies greater than that in bulk metals; while on C 60 , Ca has a binding 
energy of 1.3 eV, significantly smaller than that in bulk Ca (1.88 
eV) 18,22 , therefore metal segregation is a severe problem during 3D 
assembly. 

iii) MCOFs have intrinsic well-ordered pores freely available for 
H 2 diffusion during charge/discharge processes; while Ca 32 C 60 
material does not have this advantage. 

The porous MCOF-1 structures with stable, active metal sites are 
ideal candidates for H 2 storage. We employ first-principles DFT to 
study the interaction between MCOF and H 2 molecules, to eliminate 
uncertainties in H 2 binding configuration and energies in empirical 
calculations. Structure optimization with a full coverage of four intact 
H 2 per metal atom was performed for CaCOF-1 and LiCOF-1. We 
found that LiCOF-1 fails to bind H 2 strongly. The H 2 binding energy 
(E h ) varies from ~ 0.05 eV/H 2 to 0.02 eV/H 2 for a single H 2 and for 4 
H 2 adsorption per Li, due to the fact that Li is strongly bound to the 
upper and lower COF-1 layers and the interlayer spacing is rather 
limited (3.7 A). On the other hand, CaCOF-1 can serve as an 
excellent H 2 storage material with the H 2 binding energy of 0.13 
eV/H 2 with generalized gradient approximation (GGA) exchange- 
correlation functional and 0.42 eV/H 2 with local density approxi- 
mation (LDA), corresponding to a gravimetric storage capacity of 
4.54 wt% and volumetric density of 35.2 g/L at 0 K, very close to the 
2017 DOE target. Using VDW-DF the H 2 binding energy to CaCOF- 
1 is 0.18 eV/H 2 , a reasonable value lying between GGA and LDA 
numbers, which are known to underestimate and overestimate E h , 
respectively 11 . Simpler Grimme's D FT + dispersion approach also 
gives reasonable values for VDW interaction energies between H 2 
and metal sites in metal organic frameworks 41 . In high-level MP2 and 
quantum Monte Carlo calculations a binding energy of ~ 0.2 eV was 
obtained for H 2 binding to partially charged Ca, with basis-sets 
extrapolated to complete basis limit 42,43 . This number is close to 
our results. Based on the calculated H 2 binding energy, it is shown 
later the high H density persists at near-ambient conditions. For 
comparison, H 2 storage capacity of the COF-1 without metal inter- 
calation was found only 1.28 wt% at 77 K and ambient pressure, and 
0.08 wt% at 298 K and 20 bar in experiment 32 . 
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The nearly sixty times increase in storage density under near 
ambient conditions in COF-1 upon Ca intercalation is attributed 
to the presence of dense, active metal sites and the ideal binding 
energy to H 2 . In CaCOF-1 framework, the empty Ca 3d orbitals 
interact strongly with thep orbitals of COF-1, presenting large den- 
sity of states (DOS) around the Fermi level, see Fig. 2. These states are 
saturated upon H 2 binding. The energy distribution of H 2 orbitals is 
highly overlapped with Ca 3d orbitals, suggesting H 2 binding onto 
CaCOF-1 via a Kubas mechanism, namely, hybridization between 
the antibonding o* orbital of H 2 and the d-orbital of Ca 3,22 . 

Optimizing MCOFs: the organic building blocks. Although 
CaCOF-1 has shown outstanding property and great promise for 
H 2 storage, it is unclear whether this approach can be extended to 
other systems, or can be further improved to produce new materials 
satisfying DOE target of a 5.5 wt% density. We then explore a range 
of possible alternatives of organic building units and illuminate the 
systematic trends in H 2 binding behaviors, in the hope to optimize 
MCOF building blocks for enhanced H 2 storage capacity and to 
identify important parameters for MCOF design. 

Benzene as the reference organic unit. To seek an optimal organic 
building block, the organic unit of COF-1, a benzene ring, was first 
studied as a standard reference. Ca binding to C 6 H 6 produces a small 
energy, 0.07 eV/Ca, less than the binding energy of Ca to the single 
layer COF-1, due to limited electron conjugation in isolated benzene 
molecules. The first H 2 adsorbs onto CaC 6 H 6 in intact form, with H- 
H bond length (d HH ) expanding from 0.749 to 0.785 A. With the 
increase of the number of adsorbed H 2 molecules (n), H 2 binding 
energy and H-H bond length both first increase, then decrease, 
shown in Fig. 3. The maximum averaged H 2 binding energy E h = 
0.34 eV/H 2 is reached when n = 4; and the maximum d HH reaches 
0.804 A at n = 3. The maximum is not for n = 1 because of the 
competition between electron polarization and Pauli repulsion. We 
note that there is a strong correlation between average E h and d H H> 
which might be used to steer optimal design of H 2 storage material. 
PDOS analysis for 1-5 H 2 adsorption onto CaC 6 H 6 clearly indicates 
there is a largest mixing between H 2 states and Ca d- orbitals and a 
high occupation of H 2 antibonding states and Ca d-orbitals for n = 3 
and 4. The sixth H 2 , initially placed close to Ca, is repelled to a 



distance larger than 6 A, suggesting the system with five H 2 to be 
saturated. The maximum of five H 2 molecules offer 10 electrons, 
together with 2 electrons from Ca and 6 7i-electrons from benzene, 
the 18-electron rule is nicely satisfied as proposed for H 2 storage on 
transition metal complexes 15 . The average binding energy is ~ 0.29 
eV/H 2 . For comparison, on Ca-decorated terephthalic acid, H 2 
binding energy was calculated to be 0.20 eV in GGA and 0.24 eV 
in MP2 with completed basis set 42 . Among the adsorbed H 2 mole- 
cules, we observe that H 2 prefers to adsorb in the neighborhood of 
Ca-C 6 H 6 , than on the top site of Ca. This is also indicated by the 
difference in binding energies (0.11 vs 0.06 eV) and H-H lengths 
(0.785 vs 0.752 A) on the two sites, respectively. It is because the 
neighboring four molecules interact with Ca via the Kubas 
mechanism, as shown in Fig. 4; while the top H 2 is bound via the 
long-ranged electrostatic attraction since no overlap between H 2 
bound state and Ca d-orbital was observed in PDOS 12 . Notice that 
in the case of Ca coated fullerene 22 , the maximal number of adsorbed 
H 2 molecules per Ca is five with the same trend for E h . The Ca 
binding energy increases to 0.73 eV/Ca with two Ca bound on 
both sides of benzene, perhaps due to the more polarized n orbitals 
of benzene than the case of single Ca adsorption. Detailed analysis 
shows that upon two Ca adsorption, there are +0.99 electrons 
transferred from each Ca to benzene, larger than that for a single 
Ca adsorption (+0.88 electrons). The H 2 can be adsorbed in total, 
leading to H storage capacity of 1 1.2 wt% (0 K), with d Hli = 0.785 A; 
the E h of 0.23 eV/H 2 falls in the ideal range of 0.2-0.6 eV/H 2 . 

Comparison to other small molecules. Then a series of Ca- 
decorated small organic molecules including C 2 H m (m = 2,4,6) 
and ring structures are investigated. Although some molecules are 
not yet directly used as COF building blocks, exploring their 
potential performance yields valuable information for comparison 
to currently available materials and provides new insights for future 
applications. Compared with C 6 H 6 , we find Ca binding to the first 
three molecules is very weak, with E Ca < 0.01 eV, making them 
unstable for hosting H 2 . The binding energy of Ca to neutral C 5 H 5 
is 3.24 eV, showing a strong interaction between them. This is a result 
of electron deficiency in C 5 H 5 : to form a stable aromatic system 44 , 6 n 
electrons are required; where C 5 H 5 supplies only 5, with another 
electron from Ca. Therefore the CaC 5 H 5 complex formed is a 
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Figure 3 | (a) Average H 2 binding energy (E b ) and H-H bond length (d HH ) as a function of the number of adsorbed H 2 on Ca-C 6 H 6 . (b) Optimized 
geometry of ten adsorbed H 2 on Ca 2 C 6 H 6 . 
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Figure 4 | (a) Side and (b) top view of charge density difference after Ca binding onto benzene. Isosurfaces with values of ± 0.005 e/A 3 are shown. 
Red and blue clouds correspond to electron depletion and accumulation, respectively, (c) Same as (b) but for both two Ca and 10 H 2 binding to benzene. 
(d),(e) Same as (a) and (b) but for Ca binding on five-membered C 5 H 5 . (f) The H 2 binding energy as a function of distance (z) between Ca and H 2 on C 5 H 5 . 



product of reduction -oxidation chemical reaction, namely, C 5 H 5 is 
reduced by an electron from Ca to form a stable compound. Very 
surprisingly, the Ca-C 5 H 5 seems to be saturated and does not bind 
any H 2 . No matter where the H 2 molecule is initially placed, it is 
repelled to a distance > 4 A from Ca with negligible binding energies 
(< 0.005 eV/H 2 ), shown in Fig. 4(f). It implies the electrostatic force 
is dominant. This behavior is in contrast to that for Ca-C 6 H 6 , 
suggesting that strong interaction between Ca and the organic 
component does not necessarily introduce a strong binding to H 2 ; 
instead, the metal -organic complexes get saturated and become inert 
to H 2 binding. Therefore organic systems with too strong or too weak 
Ca binding fail to bind H 2 well. A good balance between stability (Ca- 
binding) and reactivity (Ca-H 2 interaction) is required for optimal 
H 2 storage. 

Next, we substitute carbon atoms in benzene by boron or nitrogen, 
to check whether H storage properties are affected. Typical structures 
considered and results are shown in Fig. 5 (first row) and Table 2. We 
found substitutional structures with more than three B deform and 
are unstable, thus are not listed. Ca-doped C 4 H 6 B 2 and C 4 H 4 N 2 show 
much enhancement for the first H 2 adsorption, so we study multiple 
H 2 adsorption on them. The two can both host up to five H 2 . The 
corresponding E h = 0.22 eV and d HH = 0.771 A for H 2 on 
CaC 4 H 6 B 2 , and E b = 0.14 eV and d HH = 0.784 A for CaC 4 H 4 N 2 . 
While they retain similar H 2 binding properties as that for Ca-C 6 H 6 , 
the Ca-binding is significantly improved, leading to more stabilized 
material and benign thermodynamical properties without metal 
clustering. This is due to the fact that the Ca-C 4 H 6 B 2 complex is 



formed by chemical redox reaction, where Ca electrons are trans- 
ferred to two-electron deficient C 4 H 6 B 2 molecule. If we dope two Ca 
to C 4 H 6 B 2 on both sides, ten H 2 are adsorbed, with £ Ca = 2.56 eV/Ca, 
4h — 0.771 A, and E h = 0.17 eV/H 2 , corresponding to H 2 storage 
density of 11.5 wt% (0 K). The experimental route for B-substitution 




Naphthalene 



Anthracene 




p-Terphenyl 



Diphenylethyne 



Figure 5 | A variety of organic molecules serving as COF building units 
studied in this work. Cross symbols indicate stable positions for Ca 
binding. 
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Table 2 | Ca binding energy (E Ca ), adsorbed H-H bond length (c/hh) and energy (EJ for H 2 on six- and five- membered rings 

C^H^ C4H4N2 C5H6B C4H6B2 C3H6B3 C4H4O C4H4S C4H5N 



E Ca (eV) 0.07 0.73 3.02 3.13 4.80 0.19 0.11 0.25 

E b (eV) 0.11 0.15 0.01 0.35 0.02 0.20 0.27 0.35 

^hh(A) 0.785 0.774 0.753 0.805 0.755 3.778 3.846 3.680 



might be complicated and their stability requires further investiga- 
tion. Similar results 36 ' 45 ' 46 were observed for Ca doping in B- 
substituted metal-organic framework MOF-5 and COF. 

We also consider substitution of C in the five membered rings. 
Three existing heterocycles, furan (C4H4O), tetrahydrothiophene 
(C4H4S) and pyrrole (C 4 H 5 N) are tested. Ca binding energies are 
0.10-0.25 eV on these structures, and the first H 2 is dissociatively 
adsorbed on the top of Ca, with d HH > 3.680 A (Table 2). Less than 
two more H 2 molecules can be adsorbed further, with the binding 
energy of 0.1-0.2 eV/H 2 . They are not suitable for H 2 storage. 

Large polycyclic aromatic hydrocarbons. Polycyclic aromatic 
hydrocarbons are formed when two or more benzene rings are 
connected in the same aromatic plane. The interaction energy 
between Ca and the isolated compounds increases to around 1.1 
eV, shown in Table 3. In this case, the distance between adjacent 
Ca atoms, d Ca -cc» is an important parameter. If the distance is too 
small, it could lead to metal clustering or H 2 bound to two Ca atoms, 
both reducing H storage capacity. We found on p-terphenyl (C 18 H 14 ) 
^ca-ca — 4.105 A, and four H 2 molecules between adjacent Ca 
dissociate by the strong reaction to Ca. With the distance extends 
to 6.847 A on diphenylethyne, dissociation of H 2 molecules dis- 
appears and all H 2 are kept intact with ^hh — 0.782 A. As a result, 
the total number of twenty H 2 are adsorbed on Ca-doped dipheny- 
lethyne (five per Ca), with a capacity of 10.6 wt% (0 K). However, for 
these isolated molecules, Ca binding energy is smaller than its 
cohesive energy in bulk, the metal clustering problem might occur 
if using these isolated compounds for H 2 storage. 

Optimizing MCOFs: effects of the linker groups. Besides the 
organic units, linker groups play an essential role in forming COFs 
via connecting organic units studied above. Two common covalent 
linkers, -COOH and -B(OH) 2 , are usually employed. The carboxylic 
acid groups form hydrogen-bonded assemblies 28 , while -B(OH) 2 
groups produce a six-membered ring (B 3 0 3 ) after sublimation of 
three water molecules 27 . Five-membered B0 2 C 2 ring could also 
form if hydroxyl groups (-OH) are readily available during COF 
synthesis 27 . 

To investigate the effect of linking groups on H storage, we add 
-COOH and -B(OH) 2 to four Ca-decorated compounds, including 
benzene, biphenyl, p-terphenyl and diphenylethyne. Similar to Ca- 
coated benzene, ten H 2 molecules can be stably adsorbed on 
Ca 2 C 6 H 4 (B(OH) 2 ) 2 . Slight changes in binding energy of H 2 from 
0.23 to 0.21 eV, and in d HH from 0.785 to 0.787 A are observed. Ca 
binding is also improved by 0.203 eV, due to n electron derealization 
from phenyl rings to more electron-negative linker groups. For Ca- 
coated biphenyl, we compare H 2 adsorption on compounds with 



-B(OH) 2 and -COOH linkers, shown in Table 4. While the former 
shows little difference from C 12 H 10 in every respect, the latter 
enlarges Ca-Ca distance from 3.658 to 4.155 A, leading to the number 
of adsorbed H 2 increasing from 12 to 16. This is attributed to a 
stronger interaction between Ca and phenyl rings when -COOH 
linkers are connected to phenyl rings. Therefore, the presence of 
-B(OH) 2 linkers preserves H storage properties of the original com- 
pound, while the -COOH groups increase the H 2 capacity by elong- 
ating Ca-Ca distances. 

We adopt this principle to select linker groups for Ca-decorated 
PAHs. The Ca-Ca separation on C 18 H 10 (COOH) 4 is only a little 
larger than that on C 18 H 14 (from 4.105 to 4.244 A), but H 2 capacity 
increases greatly, from 18 to 24 H 2 molecules per unit, shown in 
Fig. 6(a,b). On the other hand, since the Ca-Ca separation is already 
too large on diphenylethyne (6.8 A), we use -B(OH) 2 linkers to 
preserve its H 2 storage property. As expected, all parameters such 
as £ Ca = 0.87 eV, E h = 0.17 eV, d Ca -ca = 6.9 A, and d HH = 0.783 A 
are very close to the original complex. The same number of H 2 
molecules is adsorbed. Importantly, both systems have a desirable 
H 2 binding energy of ~ 0.17 eV/H 2 . So the Ca 4 C 14 H 8 (B(OH) 2 ) 2 unit 
could potentially serve as an optimal building block promising for 
constructing MCOF to store H 2 . 

Properties of three-dimensional CaCOF-d. The optimization 
procedures described above identify Ca intercalated COF structure 
with diphenylethyne organic unit and B 3 0 3 linkers, named as COF-d 
hereafter, as most promising for H 2 storage. Similar to COF-1, COF- 
d also forms 3D MCOF structures upon intercalation with metal 
elements such as Ca and Li, which preserves the ideal H 2 storage 
property in low- dimension materials. We then calculate the 
equilibrium structure and energy storage properties of the 3D 
intercalated CaCOF-d material in a rhombic unit cell with in-plane 
dimensions of 27.24 A X 27.24 A and vary its interlayer distance. The 
unit cell contains three branches of diphenylethyne connected to the 
center boroxine ring, namely, it consists of 42 carbon atoms, 24 
hydrogen atoms, 6 oxygen atoms, 6 boron atoms and 6 calcium 
atoms before H 2 adsorption (see Fig. 7). The optimized interlayer 
spacing is found to be 4.4 A for the minimum energy configuration, 
adopting an eclipsed stacking manner to maximize interlayer 
interactions. We note the eclipsed stacking is chosen for all COFs 
with relatively large organic building blocks, with or without metal 
intercalation 25 . The Ca binding energy in CaCOF-d, with respect to 
isolated COF-d layer and free Ca atoms as a reference, is calculated to 
be 2.79 eV/Ca (see Table 1), largely exceeding the cohesive energy of 
1.88 eV for bulk Ca. This indicates Ca atoms in this assembly cannot 
form metal clusters, similar to the case of CaCOF-1. Compared 
to CaCOF-1, the Ca binding energy is only slightly decreased by 



Table 3 | Number of adsorbed Ca atoms (Nq 0 ), average binding energies per Ca (E Ca ) and distances between adjacent Ca (c/c a -Ca) ° n 
representative PAHs. n is the number of adsorbed intact H2 molecules 

Naphthalene Anthracene Biphenyl p-Terphenyl Diphenylethyne 



CiqHs C14H10 Q2H10 QsHu Q4H10 



N Ca 4 4 4 6 4 

Ec a (eV) 1.14 1.06 1.16 1.13 0.80 

cfca-Ca(A) 3.326 3.859 3.658 4.105 6.847 

n 8 11 12 18 20 
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Table 4 | Average binding energy per Ca (E Ca ) and the distance 
between adjacent Ca (c/c a -Ca) ° n biphenyl with or without linker 
groups. Also shown are total number of adsorbed H2 molecules 
(n), and their bond lengths (c/hh) and binding energies per H 2 (Eb) 





Q2H1O 


Ci 2 H 8 (B(OH) 2 ) 2 


C 12 H 6 (COOH) 4 


Eca(eV). 


1.16 


1.25 


1.46 


Cfca-Ca ( A ) 


3.658 


3.695 


4.155 




12 


12 


16 


Eb (eV) 


0.13 


0.15 


0.16 


d»H (A) 


0.765 


0.765 


0.776 



0.04 eV/Ca. Interestingly, Li binding energy in LiCOF-d is also 
0.06 eV/Li smaller than that in LiCOF-1, indicating the dipheny- 
lethyne unit is slightly less reactive than single phenyl groups 
during MCOF formation. Both CaCOF-d and LiCOF-d intercala- 
tion are stable 3D materials. 

Three-dimensional stacking in CaCOF-d does not reduce much 
H 2 storage capacity in low- dimensioned clusters and overlayers. In 
CaCOF-d, a single Ca can adsorb four intact H 2 , all parallel to the 




Figure 6 | Optimized geometry of (a) 22 H 2 on Ca 6 C 18 H 14 , (b) 24 H 2 on 
Ca 6 C 18 H 10 (COOH) 4 , (c) 20 H 2 on Ca 4 C 14 H 8 (B(OH) 2 ) 2 . 




Figure 7 | (a) H 2 storage on two-dimensional Ca-decorated COF-d layer. 
There are five H 2 molecules per Ca adsorbed in 2D CaCOF-d. The 
hexagonal unit is blown up on the right panel for clarity, (b) Three- 
dimensional CaCOF-d structure in two views. There are four H 2 molecules 
per Ca adsorbed in 3D CaCOF-d. For clarity H 2 molecules are omitted in 
the left panel. 

plane of the COF-d layer, as shown in Fig. 7b, with the average 
binding energy — 0.15 eV/H 2 in GGA form of exchange- correlation 
energy and 0.42 eV/H 2 in LDA. We believe the exact binding energy 




50 100 150 200 

Pressure (bar) 




I ^-Without Ca 

0 ■ 1 l l ■ 

20 40 60 80 100 

Pressure (bar) 

Figure 8 | (a) H 2 storage density of 3D CaCOF-d as a function of pressure 
at 77 K, 298 K, and 400 K. H 2 binding energy around 0.29 eV is used, 
(b) Comparison of storage density in CaCOF-d (excluding free-volume 
storage) and bare COF-1 (data from experiments in Ref.27) at 77 K. 
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lies between the two values, close to the averaged value at 0.29 eV/H 2 , 
based on an empirical guideline that the binding energy averaged 
from LDA and GGA values is close to the precise Quantum Monte 
Carlo result for H 2 on B-doped fullerene 11 . A storage capacity of 4.94 
wt% at 0 K is calculated. The distance between two adjacent COF-d 
layers is 5.0 A upon H 2 binding, enlarged by 0.6 A compared to H 2 
free materials. The estimation of H 2 density excludes the free volume 
storage of H 2 in the pores of CaCOF-d, which has a diameter of ~ 20 
A. If all possible adsorption sites are considered, this porous 3D 
network satisfies the 2017 DOE target for H 2 storage. 

To unveil its practical advantage as a new class of ideal H storage 
material, Figure 8 shows the storage capacity of CaCOF-d as a func- 
tion of pressure and temperature at near ambient conditions. The H 2 
capacity at a given temperature T and pressure P is evaluated via the 
relation /(occupation number of adsorbed H 2 molecules) = kTdZI 
dfi, where k is Boltzmann constant, Z is the grand partition function, 
and \x is the H 2 chemical potential 47,48 . The grand partition function 
follows, 

Z =^2Tlexp{-g l (s l - f i)/kT}, 

gi 1 

where s t is the adsorption energy per H 2 molecule when the number 
of adsorbed molecules is /, and g t is the multiplicity (degeneracy) of 
the configuration for a given /. The chemical potential of H 2 , fi, is 
given by 

l i{P,T) = l i\T)-kTln^, 

where ff(T) is standard chemical potential and P° is standard state 
pressure. Here, /n and s t are negative values. The gravimetric density p 
of hydrogen storage is then obtained as, 

P( P > T ) = ) X Pmax- 

''max 

We adopt ambient pressure n values in the literature 48 . At 298 K 
and 20 bar, we obtain 4.87 wt% storage capacity (excluding H 2 stor- 
age in the free volume), which is substantially large than experi- 
mental values for un-doped COF-1 (0.1 wt%) 32 . The storage 
capacity is higher if H 2 adsorption onto B 3 0 3 rings and in free 
volume is considered. Most of the stored H 2 (3.2 wt%) will be released 
at 298 K and 1 bar (by solely decreasing the pressure), or at 400 K and 
20 bar (by solely increasing the temperature), where about 1.7 wt% 
H 2 retains. This favors near- ambient condition operation for H 2 
storage and release, with impressive accessible H 2 uptake. Without 
Ca, the storage capacity is very low (see Fig. 8b), and near-ambient 
condition operation is not possible due to too small H 2 binding 
energies (~ 40 meV/H 2 ). The bulk MCOFs formed by Ca intercala- 
tion have advantages of being intrinsically three-dimensional, ther- 
modynamically stable, and with well-ordered and adjustable pores, 
in contrast to previously proposed Ca 32 C 60 and other Ca- doped 
nanomaterials 22 " 24 . 

Discussion 

The 3D network based on Ca-decorated carbyne was recently pro- 
posed as stable H 2 storage media, where Ca binding energy ranges 
from 1.4 to 2.0 eV/Ca depending on Ca-Ca separation 18 . The £ Ca 's 
are close to bulk Ca cohesive energy. However, the base material, the 
3D carbyne network has yet to be synthesized in experiment. Its 
existence and room temperature stability (in particular chemical 
reactivity) remain questionable. On the other hand, crystalline 3D 
networks formed by stacking COFs via van der Waals forces have 
been explicitly demonstrated 25 ; we expect Ca-decorated COF-1 to be 
readily fabricated in experiment, by Ca intercalation between weakly 
bound COF-1 layers at high temperatures (600 K) in Ar gases. Based 
on existing materials (graphite and metals), a similar approach has 
resulted the discovery of CaC 6 , a superconducting material with 



highest transition temperature of 11.5 K among metal-intercalated 
carbon materials 49 . Hydrogen storage on metal intercalated graphite 
was also explored recently 50 " 53 . Once formed, the new Ca- 
intercalated 3D COF material is thermodynamically stable (charac- 
terized by large Ca binding energy), and would be quite inert, except 
for Ca active sites. The metal sites in MCOF might be sensitive to 
water vapor contamination, as is the case for all other H 2 storage 
materials with active sites; in real applications water concentration 
must be carefully controlled during H 2 charging and discharging 
processes. Overall, these novel properties would favor MCOF syn- 
thesis and H 2 storage applications. 

In conclusion, we have designed a new class of 3D networks based 
on first-principles: metal-intercalated COFs, which have well- 
defined crystal structures and are chemically active for molecular 
adsorption. Metal atoms prefer to bind to aromatic organic units 
rather than B 3 0 3 connecting rings in MCOFs. Li, Na, K and Ca could 
be intercalated into COF-1 and COF-d layers and form thermody- 
namically stable materials, with metal binding stronger in MCOF 
than in bulk metals. For H 2 storage, Ca is the most suitable metal 
to be intercalated in the stacking structures. Important parameters, 
such as Ca-binding and separation, and trends on 5,6-membered 
rings and larger PAHs are identified. The covalent linker -B(OH) 2 
preserves nicely H 2 storage properties, while -COOH linkers adjust 
Ca separation and in turn H 2 store capacity. Systematic optimization 
of organic building units and linker groups identifies CaCOF-d as 
most promising for near-ambient H 2 storage. The porous 3D mater- 
ial, CaCOF-d, reaches ~ 5 wt% H 2 storage capacity under near- 
ambient conditions without metal clustering problems. Routes for 
experimental synthesis are discussed. We envisage that this new class 
of MCOF porous materials, with a high density of active metal sites, 
would prove useful in a variety of applications including gas adsorp- 
tion, sensors, and air cleaning. 

Methods 

Our first-principles calculations were performed within the framework of density 
functional theory (DFT) using Vienna ab initio simulation package (VASP) 54 . The 
ultrasoft pseudopotentials 55 and general gradient approximation (GGA) in Perdew- 
Wang form for exchange -correlation energy 56 are used. For the cases of hydrogen 
adsorption in MCOFs, we also use local density approximation (LDA) for compar- 
ison 57 . Although GGA tends to underestimate H 2 binding energies and LDA over- 
estimates them, the trends for comparing relative energy differences are reliable 42 ' 58 . 
Recently developed van der Waals density functional (VDW-DF) by Langreth et al. 39 
and parameterized by Klimes et al. 40 is also employed to check the energies in some 
critical cases. The binding energies and lengths reported throughout the manuscript 
are based on GGA exchange -correlation functional unless otherwise specified 
explicitly. Spin polarization is invoked whenever necessary. We found local magnetic 
moments are zero for CaCOF systems and very small (< 0.3^) for Li- COFs and 
variations in total energy are negligible. The supercells usually contain a vacuum layer 
of > 10 A and have cell dimensions of 15 A X 15 A X 15 A and 30 A X 15 A X 20 A 
for small and large finite systems, respectively. A plane wave cutoff of 400 eV and 
single Gamma for /c-point sampling are employed for isolated systems. For periodic 
systems, a rhombic 15.20 A X 15.20 A unit cell for COF-1 and 27.24 A X 27.24 A for 
COF-d is used within the COF plane, and the interlayer distance is carefully opti- 
mized. K-point mesh of 1 X 1 X 5 is used for calculating metal and H 2 binding 
properties of MCOFs. All atoms are allowed to relax until the forces on each atom 
have magnitudes less than 0.01 eV/A. 
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